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ABSTRACT 

Given a redshift survey of galaxies with measurements of apparent magnitudes, we present a 
novel method for measuring the growth rate /(0) of cosmological linear perturbations. We use 
the galaxy distribution within the survey to solve for the peculiar velocity field which depends in 
linear perturbation theory on /3 = f(il)/b, where b is the bias factor of the galaxy distribution. 
The recovered line-of-sight peculiar velocities are subtracted from the redshifts to derive the 
distances, which thus allows an estimate of the absolute magnitude of each galaxy. A constraint 
on P is then found by minimizing the spread of the estimated magnitudes from their distribution 
function. We apply the method to the all sky K = 11.25 Two-MASS Redhsift Survey (2MRS) 
and derive (3 = 0.35 ±0.1 at z ~ 0, remarkably consistent with our previous estimate from the 
velocity- velocity comparison. The method could easily be applied to subvolumes extracted from 
the SDSS survey to derive the growth rate at z ~ 0.1. Further, it should also be applicable 
to ongoing and future spectroscopic redshift surveys to trace the evolution of f(Q) to z <~ 1. 
Constraints obtained from this method are entirely independent from those obtained from the 
two-dimensional distortion of £(s) and provide an important check on f(£l), as alternative gravity 
models predict observable differences. 

Subject headings: Cosmology: dark matter - cosmology: large scale structure of the universe 



1. Introduction 

Large scale density perturbations in the Universe 
are gravitationally unstable and grow via linear theory. 
The growing mode of large scale density perturbations, 
D(a), is characterized by the more observationally rel- 
evant growth rate 



/(«) = 



dlna 



(1) 



where a is the scale factor of the Universe and is the 
matter density parameter. It is found that the growth 
index 7 = c?ln//dlnfi is very we ll approxima ted by 
7 = 0.55 + 0.05[1 + w(z = 1)] (|Linderl [20051 ) for a 
cosmological background dominated by dark energy 
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with an equation of state, P — wpc 2 . The growth 
rate is not only sensitive to the background cosmol- 
ogy, but also to the theory of gravitation invoked as 
the driver for structur e formation. Geometric R n (e.g. 
iGannouii et alj I200S ), and Dvali-Ga badadze-Porrati 
(DGP1 (|Dvali et al.l (200ol : IWeil 120081 ) gravity models 
give substantially different behaviors of /. Models 
with dark sector long range forces, in addition to grav- 
ity, even introduce a scale dependence into / (e.g. 
iKeselman et al.ll2010h . 

Here we present a new method for constraining 
/(S7) from redshift surveys of galaxies with measured 
apparent magnitudes. Redshifts of galaxies systemat- 
ically differ from the actual distances by the line-of- 
sight components of their peculiar velocities. Hence, 
the directly measurable intrinsic luminosities or abso- 
lute magnitudes of galaxies inferred from the observed 
flux using redshifts rather than distances will show 
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larger spread than the true values. 

Gravitational instability theory allows a prediction 
of the peculiar velocity field from observed galaxy dis- 
tribution given /(fi) and the biasing relation between 
galaxies and mass. The method uses this predicted ve- 
locity field to get distances for deriving true absolute 
magnitudes. Constraints on / can then be derived by 
minimizing the scatter of the estimated absolute mag- 
nitudes from a reference distribution. Since the galaxy 
distribution could be biased relative to the mass den- 
sity field, the constraints on / are degenerate with the 



assumed biasing. Adopting linear biasing, 8 g 



bS„ 



between the galaxy number density contrast, 8 g , and 
the mass density contrast, 8 m , yields constraints on 
p = f/b that are independent from those obtained 
from the a pparent anis otropy in the observed galaxy 
clustering (|Kaiserlll988l '). 

In §2 we describe the method in detail, presenting 
general expressions and deriving the relevant approx- 
imations. In §3 we offer general analytic assessments 
of the method. In §4 we apply the method to the 
2MRS of galaxies l imited to magnitude K — 11.25 
jHuchra et alj|2005h . We conclude in §5 with a gen- 
eral discussion of the results and of the prospects for 
the application of the method to future data. 

2. The method 



Th is section is largely based on iNusser et al.l 
(|201ll ). We are given a flux limited survey of galaxies 
with observed apparent magnitudes m < mi, angular 
positions, and redshifts cz (in kms" 1 ). Let r (also 
in kms -1 ) be the luminosity distance to a galaxy in 
the sample. For simplicity of notation and description 
we assume here that the distance and spatial extent 
of the survey are small so that r is well approximate 
by the physical distance. Therefore, cz = r + v where 
v — r ■ u is the line-of-sight component of three dimen- 
sional peculiar velocity u of the galaxy. The results 
can readily be extended to the general case once we 
specify the underlying cosmological model. 

The method relies on a prediction of v from the 
observed distribution of galaxies in the survey. Since 
t — cz — v this prediction allows an estimate of the 
true absolute magnitude, 

M = m - 15 - 51ogr = M - 5 log(l - v/cz) (2) 

where the measurable absolute magnitude Mo = 
m — 15 — 5 log cz is determined from observations. Be- 
cause the peculiar velocity of a galaxy is uncorrelated 
with its true absolute magnitude, constraints on the 
underlying velocity field can be derived by demanding 
that the distribution of the magnitudes, M, is consis- 
tent with a reference distribution function (i.e. lumi- 



nosity function). The equations of gravitational insta- 
bility theory relate the underlying mass density con- 
trast, 5 m , to the peculiar velocity field u. Further, we 
will use the linearized equations in which the relation 
solely depends on /(fi). For linear biasing, 8 g = b5 m , 
between mass and galaxies the appearance of / is re- 
placed by the single parameter ft = f/b. Therefore, 
the method presented here will focus on constraints on 
/3 only. More sophisticated models for the velocity field 
involving additional parameters will not be discussed 
here. In principle, the constraint can be obtained with- 
out resorting to the luminosity function simply by min- 
imizing the variance of Mo — 51og(l — v(j3)/cz) with 
respect to j3. However, much tighter constraints are 
obtained from the full distribution. We define the lu- 
minosity function, $(M), expressed in terms of the 
absolute magnitudes, as the number density of galax- 
ies per unit magnitude. The probability, P(Mo\cz,v), 
of observing a galaxy having an observed magnitude 
Mo in a flux limited sample, depends on its redshift, 
cz, and its radi al peculiar velocity , v, and is well ap- 
proximated as (|Nusser et al.|[2oTlf ) 



P(M \cz,v) 



$(M) 



(3) 



where Mi = M i — 5 log(l — v/cz) and Moi = mi — 15 — 
5 log cz. The expression is valid as long as the relative 
errors on the measured redshift are small {<j C z/cz <C 1) 
The probability distribution of the whole sample of 
galaxies is the product of the single probabilities 



P s =HiP(Moi\cZi,Vi) 



(4) 



where i runs over all galaxies of the sample. Given a 
form for <I>(M), the parameter j3 is then constrained by 
maximizing P s in which the dependence on /3 is via the 
v(/3) as inferred from the spatial distribution of galax- 
ies. In princ iple one could use the "nonpar amatric" 
fit m ethods (|Efstathiou et all 1 19881 ; iDavis fc Huchral 
I1982T ) to approximate <5>(M). However, here we will 
only apply the method to the 2MRS sample which 
is reasonably approximated b y a Schechter luminos- 
ity function |Westoverl 120071 ). Therefore, we as- 
sume <j>(M) is w ell approximated by a Schechter form 
(|Schechterlll980h 

$(M) = 0.41n(10)<E»*10 a4(Q+1)(M *- M) 

x exp(-10°' 4(M *- M) ) . (5) 

The normalization $« does not concern us here. The 
shape parameters M, and a generally depend on the 
galaxies' type, redshift and band of observation. In 
terms of the luminosity (M = —2.5 log L + const), this 
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function acquires the simpler form 

$(M(L)) = 0.41n(10)$* (j^ exp 
Inserting all this expression into (0 gives 



P(Mo\cz;v) 



0.41n(10) 



l \ 1+a -L/L, 



T{l + a,Li/L.) 



(6) 



(7) 



where L/L, = (1 - v/czfl0 OA(M *- M °\ Li/L, = (1- 

«/c*) 2 io.°- 4 < M -- w °'>. 

To summarize, given the observed absolute magni- 
tudes Mo for a sample of objects, a Schechter model for 
the luminosity function and linear theory prediction 
for the underlying velocity field v(/3), we determine /3, 
a and M* by minimizing Eq. [4] 

3. General assessments of the method 

There are three sources of error which affect the 
derivation of /3 in the method presented here, a) "shot- 
noise" error resulting from the finite number galaxies, 
b) cosmic variance due to variations in the large scale 
structure in volumes of the universe comparable in size 
to the volume probed by the redshift survey at hand 
and c) inaccuracies in the peculiar velocity reconstruc- 
tion. In this section we offer a general assessment of 
the applicability of the method to data at higher red- 
shifts. To compensate for the degrading of the signal 
(~ v/cz) with redshift, surveys with large number of 
galaxies need to be invoked. We will only consider 
here shot-noise error. Peculiar velocity reconstruction 
errors and cosmic variance (of surveys with similar vol- 
umes) scale with redshift in a similar manner to the 
signal and hence their corresponding relative error will 
not depend on redshift. 

To the list of errors above we do not add biases 
introduced by adopting a specific form for the lumi- 
nosity function since, as we will argue in £13 .21 these 
biases are expected to be insignificant. Further, the 
issue becomes completely irrelevant for future surveys 
where the number of galaxies is large enough to allow 
the use of "nonparametric fit" techniques for modeling 
the luminosity function, relaxing the need for specific 
parametric forms. 

3.1. Sensitivity to "shot-noise" and red- 
shift 

Assume a cubicQ region at high redshift cz contain- 
ing a flux limited sample of N galaxies with redshifts 



Any survey geometry with boundary conditions allowing 
robust velocity reconstruction will suffice for our purposes. 



czi « cz. From the galaxy distribution in this re- 
gion, derive the peculiar velocity field as a function 
of /3. Consider a Schechter form for the luminosity 
function and assume cz is large enough that the lim- 
iting luminosity Li of a galaxy that could be observed 
is significantly larger than L„. In the limit Li 3> 
I/,, the F function in the expression for P(Mo\cz;v) 
in (J7)) is well approximated as T(l + a, Li/L*) = 
(Li/£„) Q exp(-Li/L*) go that, 



P(M \cz;v) = 0.4m(10)-^- ] 



(8) 



This approximation will allow an analytic expression 
for /3 by minimizing — lnP s = — ~}2\nP(Moi\czi, v(/3)) . 
For simplicity, we further approximate Vi(fi) = 
F{j3)v\i where Vi is the line of sight velocity recon- 
structed with /3 = 1. For linear velocity reconstruction 
from the galaxy distribution in real space F(/3) = /3. 
But for reconstruction from redshift space data a 
significantly bette r scaling is F = 2.5/3/(1 + 1.5/3) 
(|Davis et al.ll201lh . This is not an exact result, but it 
suffices here since we are only interested in a general 
assessment of the expected error on /3. Since /3 appears 
only via F(/3), we will perform the minimization with 
respect to F and write the result in terms of /3 at the 
end of the calculation. Using L = (1 — Fvn/czi) 2 Lo 
and Li = (1 — Fvu / 'czt) 2 L i we get 

<91nP s 



dF 
_d_ 

dF 



(9) 



E 



In 1 - F 



ALi vu 

L* CZi 



Vli 
CZ, 



F 



1 - F 



Vli 

cz, 



ALi 



where ALi = Loi—Lon and in the last step we have ne- 
glected 0(v\jcz) terms and assumed that the Hubble 
flow-like ~}2vn/cZi is negligible compared to the other 
terms. The la shot-noise error on F is, therefore, 



SF = 



d 2 \nP/dF 2 

E 



1/2 



(10) 



1 + 



ALi 
L, 



-1/2 



This expression can be easily estimated when the lu- 
minosity, Lo, is computed from the actual distances, 
i.e. czi = Ti. This means that Lo is the true intrinsic 
luminosities and, therefore, AL and vi axe uncorre- 
cted. Further, in the limit L; 3> P*, the average 
< AL/L* > is unity. From F = 2.5/3/(1 + 1.5/3) we 
get, 5/3 w 0.4(1 + 1.5/3)<5F. With all this, equation {TOj 
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gives 



5/3 = 0.4(1 + 1.5/3)(2iV)- 1/2 -^ (11) 
= 0.044(1 + 1.5/3) 



10 s \ 1/2 z 600 



N J 0.1 (Tvl 

where z =< l/z 2 >~ 1//2 ,and a v i is the rms valu e of vi 
in kms -1 . The analysis of iDavis et al.l (|201ll ) gives, 
on large scales, a v \ ~ 600 kms -1 for j3 — 1, at z — 0. 

The scaling S/3 oc N~ 1/2 z in Eg. [TT1 should be ap- 
proximately valid also when the condition Li 3> L, 
is not strictly satisfied. For 2MRS this scaling gives 
roughly the same error on /3 using distant galaxies with 
cz > 4000 kms -1 as using galaxies with the lower red- 
shifts . 

3.2. Sensitivity to assumed form of the lu- 
minosity function 

It is instructive to assess the sensitivity of /3 on 
the assumed form of the luminosity function. We 
give here a simple example in which the assumed lu- 
minosity function differs greatly from the true form, 
yet the estimate for /3 is unbiased. Consider an ideal 
volume-limited sample of galaxies with a true lumi- 
nosity function of a Schechter form with a = 0, i.e. 
an exponential distribution. Let is try to recover f3 
assuming a gaussian form for the luminosity function, 
P(M a \cz,v) oc al x {l - Fi> 1 ) 2 e- (i ~ im)2/2<T S where, 
like in the previous section, v(/3) = F{j3)v\. Minimiz- 
ing the quantity, — \nP(Moi\cZi, Fvu), with respect 
to Lm, &l and F, yields, respectively, 



— L r . 



< L > , 



(12) 



= 



o'l- < (L 
E 



o 2 >, 

Lo(Li — L„ 



Vl 



CZi CZi 



F- 



For simplicity we further assume that Lo is computed 
with CZi = Ti so that the true solution ft — (the 
generalization to cases where the true j3 is different 
from zero is trivial). Hence, Lo are equal to the true 
luminosities and, therefore, follow an exponential dis- 
tribution for which < Lo >= L» and < L 2 >= 2L 2 . 
Further, there is no correlation between Lo and — , 
meaning that the average of products of powers of 
(1 — Fvu/czi) and Lo is the product of the averages. 
In the limit N — > oo, straightforward algebraic manip- 
ulation then yields F(/3) = 0, L m = L„, and a\ = L 2 . 
Therefore, in this example, where the assumed and 
true luminosity functions differ grossly, the best fit /3 
is unbiased. This is not surprising since the underlying 
principle of the method is a reduction in the spread 



of L. Therefore, the assumed form of the luminos- 
ity function should only affect the weighting given to 
galaxies in a certain luminosity range rather, than the 
best fit p. Assuming a wrong luminosity function in- 
creases the random error on j3 but does not introduce 
any systematic biases. 

4. Application to 2MRS 

In this section we apply the method outlined above 
to the all sky 2MRS consisting of 23,200 galaxies down 
to the magnitude K — 11.25. Details about the cat- 
alog, including the precise completeness, sky coverage 
and selection effects can be found in (|Huchra et al.l 
|2005| ). The preparation of the catalog for the purpose 
of the applicat i on of the method is is done similarly 
to IDavis et al.l (|201ll ) . The peculiar velocity field is 
derived from the galaxy distribution in the 2MRS for 
an array of j3 values using the linear theory method- 
ology of iNusser fc Davis! (| 19941 ') . The derived veloc- 
ity field is robust within cz < 10 4 kms -1 , above 
that redshift discreteness effects become important. 
Hence we limit the analysis to the 18, 000 galaxies 
with cz < 10 4 kms -1 . In the derivation of the veloc- 
ity fields, the galaxy distribution is smoothed with a 
gaussian window of constant width of 400 kms -1 . To 
further remove strong nonlinearities the derived three- 
dimensional velocity fields are smoothed with an Gaus- 
sian window of constant width, R s . Linear theory re- 
covers the flow p attern reasonably w e ll even at S ^ 1, 
but n ot beyond (|Nusser et al.l Il99ll ; iBranchini et al.l 
l2002al ). Therefore, although the peculiar velocity field 
is predicted from the distribution of all galaxies, in 
the maximization of P s to assess the robustness of 
the method we remove galaxies in regions with density 
contrast higher than S cut as listed in Table [T] for both 
values of R s . Using the expression (J7J we minimize 
— lnP s = — £\ lnP(Moi\czi\ v((3)) (the summation is 
over all galaxies) with respect to j3 and the Schechter 
parameters, a and M». 

4.1. Error estimation based on mock 2MRS 
catalogs 

The overall expected errors in /3, including possible 
biases, are based on mock catalogs designed to match 
the general properties of the 2MRS. For this purpose 
we use 135 2 MRS mock catalog s very similar to those 
compiled by IDavis et al.l (|20 1 1 ) . These catalogs are 
extracted from a parent mock catalog of the Two M i- 
cron All Sky Survey (2MASS) (|Skrutskie et alj|2006j ). 
The parent catalog is generated from the Millennium 
simulation (jSpringel et al.l 120051) using semi-analyti c 
galaxy formation models (|De Lucia fc Blaizoj l2007h . 
All 135 2MRS catalogs satisfy the following conditions. 
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a) The central "observer" in each mock is selected to 
reside in a galaxy with a quiet velocity field within 
500 kms^ 1 , similar to the observed universe, b) The 
motion of the central galaxy is 500 to 700 km s" 1 , and 
c) The density in the environment of the local group, 
averaged over a sphere of 400 km/s radius , is less than 
twice the normal. These conditions select central ob- 
servers that are similar to the conditions of our own 
Local Group. 

The luminosity function in each mock is approx- 
imated by a Schechter form. Using the galaxy dis- 
tribution in each mock, the corresponding peculiar 
velocity field is generated for an array of ft values. 
The mean and standard deviation of those best values 
from the 135 catalogs is /3 = 0.49 ± 0.13. The cos- 
mological density parameters of the Millennium sim- 
ulation are fl = 0.25 and A = 0.75 for matter and 
cosmological constant. This yields / = f2 0,55 = 0.47 
(|Linderl I2005T ) . A calculation of the rms galaxy den- 
sity fluctuations of the mocks yields b = 1 ± 0.1 for 
the mocks. The la rms scatter of the best fit P val- 
ues is 0.13. This scatter reflects "shot-noise" errors 
due to the finite number of galaxies, cosmic variance 
due to the limited volume covered by the 2MRS, and 
inaccuracies in the reconstruction of the peculiar ve- 
locity field by means of linear theory. Shot-noise is 
subdominant when the method is applied to all 2MRS 
galaxies within 10 4 kms -1 . Since P s contains no in- 
formation about cosmic variance and reconstruction 
errors, the width of — 21nP s versus p, reflects shot- 
noise errors only. From the width of — 21nP s we get 
la shot-noise error of 8f3 sn m 0.055 in the application 
to 2MRS within cz < 10 4 kms -1 . Cosmic variance 
error could be estimated by running the method with 
the actual velocity field of galaxies in the mock. This 
amounts to 5/3 cv as 0.09. Adding errors in quadratures, 
we infer a velocity reconstruction error of 8p vrec fa 0.1, 
comparable to Sp cv . 

We will apply the method to various cuts of the 
2MRS. All corresponding errors listed in Table [T] are 
based on similar cuts taken from the mocks. Shot- 
noise errors, Sf3 sn , will be treated as (nearly) indepen- 
dent of P (see t|3.ip . while cosmic variance, Sp cv and 
velocity reconstruction errors, S/3 vrec , are assumed to 
be proportional to (3, because of the (3 dependence in 
the reconstructed field. 

4.2. /3 from the real 2MRS 

The results are presented in Table [1] and the top 
panel of Fig. [T] The plotted quantity, Ax 2 , is — 21nP s 
evaluated as a function of P minus its value at the 
minimum. In the figure, the parameters (q,M«) are 
fixed at their best fit values. The overall errors listed 
in Table [T] are based on the analysis of the mocks 



in £14.11 The la shot-noise, S sn , given by the width 
of the curves at Ax 2 = 1 are significantly smaller 
than the corresponding total errors quoted in the Ta- 
ble [T] i.e. 8 f} a n is not the main source of error for 
the 2MRS catalog. Table [T] lists results only for R s — 
Gh^ 1 Mpc, however, an application of the method with 
R s = 10/i -1 Mpc and 12ft" 1 Mpc yields consistent re- 
sults within the la errors in the table. T he derived P 
for bo th choices of 5 cut are consistent with lDavis et al.l 
(l201lh and the values of a and M* agree well with 
IWestoverl (|2007l ). 

It is useful to examine how the best fit P changes 
when in the maximization of P 3 we include only dis- 
tance galaxies with cz > 4000 kms -1 . The results are 
given in the Table [1] and in the bottom panel of Fig. [T] 
showing Ax 2 versus P for this case. Shot-noise as in- 
dicated by width of curves increases with respect to 
the full sample (see top panel of the same figure). The 
constraints from the distant cut are still reasonably 
tight and, within the total la errors, are fully consis- 
tent with those obtained from the whole sample. 

We also applied the method to galaxies in the 
Northern and Southern Galactic hemispheres, sepa- 
rately. The results, listed in the Table [1] show a 
non-negligible difference between the derived P in the 
two hemispheres. The corresponding error on P is 
based on the application of the method to "North- 
ern" and "Southern" hemispheres in the mocks. We 
attribute the difference between the North and the 
South to cosmic variance and confirm that it is con- 
sistent with the mocks. For each mock catalog we 
compute Pnorth and p SO uth and find an rms value 
<(/3„o rtll -& outl J 2 > 1/2 =0.2. 

5. Discussion and Conclusions 

We have presented a new method to determine P 
from galaxy redshift surveys. The method is entirely 
independent of distance indicators such as the Tully- 
Fisher relation and of anal yses of aniso tropic correla- 
tion functions, £(r p , 7r) fe.g. lKaiserll 19871 ). As a prelim- 
inary application of the method, we have resorted to 
the K=11.25 flux limited 2MRS all sky survey. In the 
maximization procedure of P(Mo\cz, v(P)), galaxies in 
dense regions should be excluded since the flow pattern 
in these regions is not well recovered by linear theory. 
For our adopted density contrast cut S cut = 1, we get 
a best fit P = 0.35 ±0.1 for velocities smoothed with a 
gaussian window of 6/i -1 Mpc in width. However, the 
results obtained with 8 ou t = 2 and 5 cu t = co are con- 
sistent with this best fit value and are reported in the 
Table [Tl These constraint s on P agree very well with 
those of lDavis et all (|201ll ) who compared the peculiar 
velocities of the S piral Field I Band (SFlH — h) catalog 
of spiral galaxies (|Masters et al.|[2006l ; ISpringob et al.l 
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Table 1: Derived parameters, (3 = f/b, a and M* for various cuts of the 2MRS. Galaxies in regions with 
density contrast larger than 6 cu t and redshifts less than cz cut are excluded from the maximization of P s . The 
quoted errors on f3 are based on mock catalogs and include shot-noise, cosmic variance and inaccuracies in 
the peculiar velocity reconstruction. The results given here correspond to R s = 6ft. -1 Mpc. 
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120071 ) with the velocity field predicted from the 2MRS. 

There are three sources of uncertainties which con- 
tribute to the error budget on the estimated /?: 

(1) "shot-noise" due to the finite number of galaxies 

(2) cosmic variance which reflects the variation of the 

large scale structure in random volumes compa- 
rable in size to the volume covered by the data 
set under consideration. 

(3) inaccuracies in the linear methodology for recon- 

structing the peculiar velocity from the galaxy 
distribution. 

For the 2MRS sample of ~ 18, 000 galaxies within 
cz — 10 4 kms -1 , shot-noise is subdominant. In- 
creasing the number of galaxies in the sample without 
probing larger volumes will not tighten the constraint 
significantly. Cosmic variance can only be reduced 
by pushing toward deeper and larger surveys. The 
main ga laxy sample of the Sloan Digital Sky Survey 
(SDSS) (jStrauss et alJl2002l ) already offers this oppor- 
tunity. Let us consider the ~ 7500 deg 2 patch around 
the Northern Galactic cap in the SDSS-DR7 release 
(jAbazaiian et al-lliou^ . A shell with Az = 0.1 cen- 
tered at z ~ 0.135, close to the peak of the galaxy 
dN/dz has a comoving volume ~ 20 times larger than 
that of the 2MRS and it could be divided into ~ 12 
independent cubes of 200/i~ 1 Mpc, each containing as 
many galaxy as 2MRS. Applying our method to each 
of them would dramatically decrease cosmic variance 
and errors in the reconstructed velocities. According 
to Eq. [11] the now dominant shot-noise error would 
be 5/3 ~ 0.06, twice as small as in the 2MRS case, as 
shown in Tabled 



Shot-noise errors increase linearly with redshift and 
they will be the limiting factor in limiting the preci- 
sion of constraining ft from the method presented in 
this paper. Next-generation, large redshift surveys like 
Boss, BigBoss and EUCLID will allow an application 
of our method to measure /3 out to z ~ 1. We have ap- 
plied Eq. [TT]to estimate the expected shot noise error 
on f3 for all these surveys. The results are listed in Ta- 
bled The ongoing SDSS-III (BOSS) survey will target 
highly luminous galaxies with a nearl y constant num- 



ber d ensity n ~ 3 • 10 4 h 3 Mpc 3 ( Eisenstein et al] 



1201 lh over the redshift range [0.2 — 0.6]. Because of the 
relatively low number density of objects our method 
will measure /3 with uncertainties twice as large as 
for the 2MRS case. The BigBOSS survey will expand 
Boss, both in sky coverage and depth. The expected 
number of ELGs (|Schlegel et al.1 12011T ) will be large 
enough to decrease the errors on f3 significantly. In 
fact, the reduction of the relative errors will be even 
larger since, in this redshift rage, j3 is an i ncreasing 
funct ion of z. Finally, the EUCLID survey (|Laureiisl 
l2009h will reduce errors even further in the redshift 
range z = [0.7 — 1.0]. 

How do the f3 estimates obtained with our method 
compare with those obtained from the analysis of the 
anisotropy pattern of the two-point correlation func- 
tion C( r P) 7r ) ? At low (2 < 0.2) reds hift the £(r p ,7r) 
meth od has be en applied to SDSS jTegmark et al.l 
l2006h and 2dF jHawkins et al.1 |2003t ) allowing to es- 
timate f3 with an error Sf3 — 0.1 — 0.15. The appli- 
cation of our method to 2MRS already gives /3 with 
a similar preci sion and the upcomi ng application to 
SDSS-II data (|Bernardi et alj l2003h will allow us to 
estimate j3 with a precision sufficient to test the va- 
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0.25 0.3 0.35 0.4 0.45 0.5 0.55 

P 

Fig. 1. — Top: Curves of A\ 2 as a function of (3 
for R s = 6h- x Mpc and 10ft," 1 Mpc, and 5 cut = 2 
and 1, as indicated in the figure. The width of 
the curves reflects uncertainties due to the finite 
number of galaxies ( "shot- noise" ) and do not in- 
clude cosmic variance and inaccuracies in the pe- 
culiar velocity reconstruction. Bottom: The same 
as the top panel, but excluding nearby galaxies 
with cz < 4000 kms -1 in the minimization of P s . 
The widths of curves here are substantially larger 
than in the top panel due to the smaller number 
of galaxies and the larger redshifts. 



lidity of popular alternative gravity model, lik e the 
5-dim ensional brane- world of iDvali et al. (2000); IWeil 
(|2008l ) as well as the possibility of a coupling be- 
tween the dark energy and the dark matter sectors 
(|di Porto fc Amendolail2008l ). 

The first measurement of the growth rate a t large r 
redshifts has been performed by iGuzzo et al.l (|2008). 
From the observed n) of VVDS galaxies they ob- 
tained p(z = 0.77) = .70 ± 0.26. Mor e recently, 
the Wiggle-z experiment jBIake et al.ll201lh has mea- 
sured the normalized growth rate /(n)crg in 4 red- 
shift bins to z — 0.9 with an error 5f(Q,)as ~ 0.1. 
This is already comparable with the expected perfor- 
mance of our method on future datasets. Indeed, our 
method does not compare favorably with the analy- 
sis of galaxy clustering in future surveys. In the case 
of EUCLID, the goal is to estimate the growth rate 
from £(r p ,n) with a precision of 0.01 at z = 1, if 
the rms mass fluctuation erg c an be determined ac- 
curately (jSong fc PercivaJ 120091 ) . In this case our al- 
ternative method to measure f3 will constitute an ef- 
fective way to keep systematic errors below 5/3 = 0.1. 
It is important to note that if galaxy bias will only be 
constrained at the 10% level then both methods, the 
analysis of £(r p ,n) and the one proposed here, will 
constrain the growth rate /(fi) with similar precision. 
Further, on-going and planned redshift surveys will 
deliver velocity dispersion for all the elliptical galax- 
i es will be obtained. Us ing the Faber- Jackson relation 
(|Faber fc Jacksor]|l976h . between luminosity and ve- 
locity dispersion of elliptical galaxies, in conjunction 
with our method will produce even tighter constraints. 
Additional constraints on / could also be obtained 
from the expected large scale supernova survey (e.g. 
iBhattacharva et al.|[201lT ). 

Quasi-linear (i.e. the mildly non-linear regime) 
dynamical reconstruction methods offer a substan- 
tial improvement in the accuracy of the predicted 
peculiar velocities. In particular, t he Fast Action 
Metho d (|Nusser fc Branchinil |2000| ; iBranchini et al.1 
l2002lj) adap t ation of Peebles' least action principle 
(|Shava et all Il995f ) is fast and easy to implement. 
This method is significantly better than linear theory 
for the reconstruction of the peculiar velocity field on 
small scales and also in dense regions. The linear the- 
ory relation between mass and velocity solely depends 
on the growth rate f(Q) = . In the quasi-linear 
regime, there is an additional explicit dependence on 
fi, raising the possibility of separate constraints on 
Q and the growth index 7. However, this does not 
seem a promising route to break the f2 — 7 degeneracy 
since the explicit O depe ndence is very weak (~ Q 2 ) 
(N usscr fc Colberdll998h . 

Another improvement could be achieved by use of 
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Table 2: Expected shot-noise errors on (3 estimated from Eq.[TT]for present and future galaxy redshift surveys. 
To minimize evolutionary effects all errors are estimated in bi ns Az = 1. Th e reference /3(z) function has been 
obtai ned by consider i ng the growth rate of a ACDM model of lLinderl (|2005l ) and a linear bias b(z) = y/Jl + z) 



as m 



Di Porto et al.l (|20 1 lh . The references to the surveys' parameters are given in the text. 



Survey 


sky 


(z) 


N t 




5p 




deg 2 










SDSS-II 


7500 


0.135 


3 • 


10 b 


0.07 


Boss 


10000 


0.25 


9 • 


10 4 


0.18 


Boss 


10000 


0.35 


3 ■ 


10 5 


0.17 


BigBoss 


14000 


0.45 


5 ■ 


10 5 


0.16 


BigBoss 


14000 


0.55 


1.3 


■ 10 6 


0.14 


BigBoss 


14000 


0.65 


2 ■ 


10 6 


0.13 


EUCLID 


15000 


0.70 


3.5 


■ 10 6 


0.11 


EUCLID 


15000 


0.80 


6.5 


• 10 6 


0.09 


EUCLID 


15000 


0.90 


7- 


10 6 


0.09 


EUCLID 


15000 


1.00 


7.5 


• 10 6 


0.10 


EUCLID 


15000 


1.10 


7- 


10 6 


0.11 



nonparametric fit techniques for modeling the galaxy 
luminosity distribution. Although the Schechter form 
is a good approximation for the 2MRS, it is likely 
less successful for larger and deeper data sets. A va- 
riety of such nonparametric fit methods exi st (e.g. 
lEfstathiou et al ] ll98Sl : iDavis fc Huchrl] Il982h could 
easily be incorporated in the method presented here. 

Angular coherent photometric mis-calibrations are 
likely to contaminate the data at some level. However, 
as long as the density field inferred from the galaxy 
distribution is not significantly affected, the method 
should yield an unbiased /3. This is because the un- 
derlying velocity field should be uncorrelated with ob- 
servational mis-calibrations. If systematic biases intro- 
duce serious spurious modes in the density field then it 
is doubtful if the data could be useful for any analysis 
of large scale clustering. 
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